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ABSTRACT 

Milli-arcsecond astrometry provided by Hipparcos and by radio observations makes it possible to 
retrace the orbits of some of the nearest runaway stars and pulsars to determine their site of origin. 
The orbits of the runaways AE Aurigae and ^ Columbae and of the eccentric binary l Orionis intersect 
each other ~2.5 Myr ago in the nascent Trapezium cluster, confirming that these runaways were formed 
in a binary-binary encounter. The path of the runaway star C Ophiuchi intersects that of the nearby 
pulsar PSR J1932+1059, ~1 Myr ago, in the young stellar group Upper Scorpius. We propose that 
this neutron star is the remnant of a supernova that occurred in a binary system which also contained 
C Oph, and deduce that the pulsar received a kick velocity of ~ 350 km in the explosion. These two 
cases provide the first specific kinematic evidence that both mechanisms proposed for the production of 
runaway stars, the dynamical ejection scenario and the binary-supernova scenario, operate in nature. 

Subject headings: Astrometry — stars: early-type — stars: individual {C, Oph, AE Aur, Col, t Ori) 
~ pulsars: individual (PSR J1932-M059) 



1. INTRODUCTION 

Runaway stars (Blaauw 1961) distinguish themselves 
from the normal population of early-type stars (spectral 
type O and B) by their large pecuHar velocities (up to 200 
km s""'^) with respect to the mean Galactic rotation and 
by their isolated locations. About 30% of the O stars and 
5-10% of the B stars are runaways. Since early-type stars 
are relatively young (a few Myr to ~50 Myr old), the dis- 
tances traveled by the runaways are relatively small (a few 
hundred pc to a few kpc) and it is therefore possible, in 
some cases, to identify the parent cluster, i.e., the stellar 
group where the runaway originated (e.g., Blaauw 1993). 

Two mechanisms for the production of runaway stars 
remain viable, the binary- supernova scenario (Zwicky 
1957; Blaauw 1961) and the dynamical ejection scenario 
(Poveda, Ruiz & Allen 1967; Gies & Bolton 1986). In the 
former the runaway receives its velocity after a supernova 
explosion in a massive close binary; after the explosion the 
binary (sometimes) dissociates and the secondary starts 
to move through space with a velocity comparable to its 
pre-explosion orbital velocity. In the latter scenario the 
runaway gains its velocity through a dynamical interac- 
tion with one or more other stars. The most efficient en- 
counter is that between two hard binaries which, in most 
cases, results in the ejection of two runaway stars and one 
eccentric binary (Hoffer 1983; Mikkola 1983). 

Which of the two mechanisms dominates the creation 
of runaway stars has been debated vigorously in the past. 
The modest accuracy and incompleteness of previous data 
sets did not allow distinguishing convincingly between the 
two mechanisms: both scenarios were consistent with the 
statistical properties of the ensemble of runaway stars. 
The availability of Hipparcos milli-arcsec (mas) astrome- 
try for the nearby stars (ESA 1997), and pulsar astrometry 
provided by VLBI and timing measurements, stimulated 
us to revisit this issue, as it allows calculation of the past 
orbits of these objects with sufficient accuracy to estab- 



lish their parent group, and the mechanism that formed 
them. A past encounter between a runaway star and a 
pulsar would provide compelling evidence for the binary- 
supernova scenario, while an encounter of two (or more) 
runaway stars would strongly suggest the dynamical ejec- 
tion scenario. In the course of this work, we could identify 
the specific formation scenario for three nearby runaways 
(the pair AE Aur & fi Col and C Oph) with near certainty. 

In this Letter we summarize the results for these ob- 
jects, as together they demonstrate that both mechanisms 
for the creation of runaways stars operate in nature. De- 
tails of the orbit integrations, the Galactic potential used, 
and the simulations performed to validate the results, as 
well as a discussion of the other nearby runaway stars and 
pulsars with accurate astrometry, are reported in Hooger- 
werf et al. (2000, HBZ). 

2. A STELLAR ENCOUNTER IN ORION 

The Orion association (Ori OBI) has three known run- 
away stars: AE Aur, /i Col, and 53 Ari (Blaauw 1961). 
The first two of these form a pair; they have almost sim- 
ilar spectral types (09.5V and O9.5V/B0), and move in 
opposite directions at 100 km s~^ each, leaving Ori OBI 
about 2.5 Myr ago. These similarities led Blaauw & Mor- 
gan (1954) to suggest that the two stars were formed in 
the same event. Gies & Bolton (1986) proposed that 
AE Aur, /i Col, and the massive eccentric binary i Ori 
(09III-I-B1III) are the result of a binary-binary encounter 
which ejected the two runaways. Both stars have normal 
rotational velocities (25 and 111 km s^^ for AE Aur and 
/i Col, respectively) and AE Aur has a normal He abun- 
dance. The He abundance of ^ Col is unknown. Thus 
neither runaway shows signs of previous mass transfer in 
a close binary system, and this most likely excludes the 
binary-supernova scenario as the origin of the two run- 
aways. Dynamical ejection is by far the favorite scenario 
for the origin of the high velocities of AE Aur and Col. 

We supplemented the Hipparcos astrometry with the 
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best available radial velocities in order to retrace the or- 
bits of AE Aur, jjL Col and t Ori. We performed 2.5 x 10^ 
orbit integrations to sample the errors in the positions and 
velocities of the objects (see HBZ). The integrations show 
that ^2.5 Myr ago the three stars were very close together: 
the distribution of minimum separations obtained from the 
orbit integrations is consistent with the stars being located 
in exactly the same position in space at the same time, 
2.5 ± 0.2 Myr ago. We therefore conclude that the two 
runaway stars and the binary t Ori must once have been 
part of the same cluster, and were all ejected following a 
dynamical interaction. 

A natural question to ask is which cluster hosted the four 
stars before their encounter? Assuming that the center-of- 
mass motion of the four stars was similar to that of the 
parent cluster, and using conservation of linear momen- 
tum, provides the position and velocity of the parent clus- 
ter ^^2.5 Myr ago. We integrated its orbit forward in time 
to the present (Figure 1). The resulting properties of the 
parent cluster, in particular its distance and position on 
the sky, agree very well with that of the Trapezium cluster 
(see Table 1 and Figure 1). 

Several other properties of the Trapezium cluster 
strengthen the conclusion that this is the most likely can- 
didate for being the parent of AE Aur, ^ Col, and b Ori. 
First, the Trapezium is a very young cluster, ~2 Myr 
(Palla & Stabler 1999); its density (>20,000 stars pc'^ 
in the center) is still high enough for stellar encounters to 
occur. At the same time, the Trapezium is old enough to 
have existed when the runaways left the cluster ~2.5 Myr 
ago. Secondly, the Trapezium shows a strong mass segre- 
gation (Zinnecker, McCaughrean & Wilking 1993; Hillen- 
brand & Hartmann 1998); this concentration of massive 
stars increases the probability for dynamical encounters 
between stars. Thirdly, the binary fraction in the Trapez- 
ium is high, 60 100% (Prosscr ct al. 1994; Weigelt et al. 
1999) . This also increases the chance of dynamical encoun- 
ters since binary- binary collisions are the most efficient. 

3. A SUPERNOVA IN UPPER SCORPIUS 

Blaauw (1952) identified the 09.5V star C Oph as a 
runaway originating in the Sco 0B2 association. The star 
could either have left the Upper Scorpius subgroup ~1 Myr 
ago or the Upper Centaurus Lupus subgroup ~3 Myr ago. 
The intrinsic properties of C Oph (observed rotational ve- 
locity Vrot sin i — 350 km s~^ and He abundance Y = 0.40) 
indicate that the star previously experienced mass trans- 
fer in a close binary system. As ( Oph is single at present 
(Gies & Bolton 1986), the binary must have dissociated 
sometime in the past. This led to the suggestion that 
C Oph is a runaway created in a binary supernova explo- 
sion. Finding a compact object formed in the same event 
would prove that this assumption is correct. 

Radio pulsars are the only neutron stars for which reli- 
able and accurate proper motion measurements are avail- 
able. Of all pulsars in the Taylor, Manchester & Lyne 
(1993) catalog, only seven are within 1 kpc, and have 
proper motions with better than 10% accuracy. Six of the 
seven cannot be related to ( Oph, based on their (three- 
dimensional) position and velocity relative to ^ Oph, their 
position and velocity relative to the Galactic plane, or their 
characteristic age P/{2P) (see HBZ for details). The pul- 



sar that remains is PSR J1932+1059 (a.k.a. B19294-10); 
it has a characteristic age of ^3 Myr, and it traversed the 
Upper Scorpius region about 1 Myr ago if its (unknown) 
radial velocity Wrad = 200 ± 50 km s~^ (cf. Figure 2). 

The pulsar currently moves away from the Galactic 
plane with a z-velocity of ~40 km s~^, and is located at 
a Galactic latitude of 6 —4 degrees. Assuming that 
most neutron stars are created in or near the Galactic 
disk, PSR J1932-I-1059 must have formed either recently 
or some 50 Myr ago, when its past orbit again crossed 
the Galactic plane (e.g., Blaauw & Ramachandran 1998). 
Taking into account its characteristic age, it is natural to 
assume that the pulsar was formed recently. The only site 
of active or recent star formation along the pulsar's path 
is Upper Scorpius, and this young stellar group therefore 
is the only likely birth-site for PSR J1932-hl059. 

We integrated the orbits of C, Oph, the pulsar, and Up- 
per Scorpius back in time in a standard Galactic potential. 
We performed 3 x 10^ such integrations to sample the error 
distibutions in position and velocity. The main uncertain- 
ties are the errors in the parallax of Q Oph (7.1 ± 0.7 mas) 
and the pulsar (5 ± 1.5 mas, Campbell 1995), and the re- 
maining range in Wrad- The distribution of the minimum 
distance between C, Oph and the pulsar is consistent with 
zero distance, i.e., with both objects being in the same 
location, 1.0 ± 0.1 Myr ago in Upper Scorpius (Figure 2). 
This is strong evidence for the binary supernova scenario. 

The kinematics of the expanding HI shell which sur- 
rounds Upper Scorpius requires a supernova explosion 1- 
2 Myr ago, and the present-day mass function of the sub- 
group suggests that originally one more massive star or 
binary (~40 Mq) must have been present (de Geus 1992). 
The simplest interpretation is that this was a close binary 
containing Oph and the progenitor of PSR J1932+1059. 

The value of P/2{P) is an uncertain age indicator, and 
the - 3 Myr for PSR J1932+1059 is consistent with the 
kinematicl age of 1 Myr derived here. The implied period 
at birth is 0.18 seconds; the current period is 0.22 seconds. 

Pulsars are expected to receive a kick velocity Ukick at 
birth (e.g., Lai 2000). Observations of the ensemble of 
pulsars suggest that the magnitude of Wkick is a few hun- 
dred km s~^ (Hartman 1997; Hansen & Phinney 1997). 
Assuming that C Oph and PSR J1932-M059 originated 
in the same binary allows an individual determination of 
Vkick- Based on the magnitude of the space velocities, and 
the angle between the orbits of C Oph and the pulsar, we 
obtain |i?kick| = 350 ± 50 km s~^. New VLBI observations 
of the proper motion and parallax of PSR J1932-I-1059 are 
being obtained by Campbell (priv. comm.), and will make 
it possible to improve this estimate further. 

4. CONCLUSIONS 

The cases described in the previous two sections provide 
the first specific evidence that both the binary-supernova 
scenario and the dynamical ejection scenario produce run- 
away stars. This result is based on a detailed analysis 
of the orbits of three runaway stars, in contrast to ear- 
lier investigations which mainly focused on the statistical 
properties of a set of runaway stars. 

Determination of the fraction of the runaway stars cre- 
ated by either mechanism will put strong constraints on 
cluster formation theories and on stellar evolution theo- 
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ries. Numerical simulations such as those performed by 
Portegies Zwart (1999) show that runaway statistics can 
put limits on, for example, the number of primordial bina- 
ries. However, the current samples of runaway stars and 
pulsars with accurate astromctry arc severely incomplete. 
The Hipparcos Catalog is complete only for y < 7.3 — 9 
mag (depending on latitude), and less than a third of the O 
and B stars in it have a measured radial velocity. Because 
of beamed radio emission, we cannot observe all pulsars, 
and not all of those that do radiate in our direction have 
been found in the Solar neighbourhood; of these, only a 
few have an accurately measured proper motion. Even so, 
the available data supports the tentative conclusion that 
both mechanisms contribute about equally to the produc- 
tion of runaway stars (see HBZ). 



The approach followed here promises many new results 
once the next generation of astromctric satellites, e.g., 
FAME and in particular GAIA, are launched and VLBI 
techniques are developed fully. These will improve the 
accuracy of the astromctry for stars and pulsars into the 
micro-arcsecond regime. With a fainter limiting magni- 
tude, many more runaways will have well-determined po- 
sitions and velocities. It will then be possible to identify 
the parent clusters of these objects, and to learn much 
about the ages and kick velocities of individual pulsars. 

It is a pleasure to thank Adriaan Blaauw, Bob Camp- 
bell, Nicolas Crctton, Ed van den Heuvcl, Rob den Hol- 
lander, Michael Ferryman and Simon Fortegies Zwart for 
useful comments and suggestions. 
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Fig. 1. — Top: The past orbits on the sky of the runaways AE Aur (dotted hne) and Col (soUd Hne) and of the binary l Ori. The starred 
symbols depict the present position of the three stars. The orbit of the parent cluster is denoted by the grey solid line. The dots denote all 
stars in the Hipparcos Catalogue brighter than V = 3.5 mag (large) and with 3.5 mag < V < 5 mag (small, only the O- and B-stars) (cf. 
figure 1 in Blaauw & Morgan 1954). The Orion constellation is indicated for reference. Bottom: The predicted position of the parent cluster 
(black contours) together with all stars in the Tycho Catalogue (ESA 1997) in the field. The size of the symbols scales with magnitude. The 
black and dark grey lines are the past orbits of i Ori and the Trapezium, respectively (see top panel). The grey contours display the IRAS 
100 micron flux map, and mainly outline the Orion Nebula. 



Table 1 



Predicted properties of the parent cluster of AE Aur 
AND ji Col, and observed properties of the Trapezium 



Property 


Parent cluster 


Trapezium 


Unit 


D 


425-450 


450-500 


pc 




(83?9, -5?2) 


(83?8, -5? 4) 






(1.7,-0.2) 


(2.7,-0.9) 


mas yr~^ 


{e,b) 


(208?9,-19?2) 


(209?0, -19?4) 


mas yr""' 




(0.9, 1.4) 


(2.0,2.0) 


Vrad 


27.6 


24 


km s~^ 
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Fig. 2. — The orbits of Oph, PSR J1932+1059, and Upper Scorpius. The present positions are denoted by a star for the runaway, a filled 
circle for the pulsar, and by a filled square for the center of the association. The top panel shows the distance vs. Galactic longitude of the 
stars. The bottom panel shows the orbits projected on the sky in Galactic coordinates. The small open circles in the bottom panel denote the 
present-day positions of the O, B, and A-type members of Upper Scorpius (dc Zceuw et al. 1999). The large circle denotes the position of the 
association at the time of the supernova explosion (10 pc radius). This figure assumes a set of space motions consistent with all observables. 



